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■ Abstract. We discuss two recent topics on heavy-ion quasi-elastic scattering at energies around 

I the Coulomb barrier. The first topic is an application of quasi-elastic scattering at deep-subbarrier 

. energies to extracting the surface diffuseness parameter of the nucleus-nucleus potential. The sec- 

^ ■ ond topic is a coupled-channels analysis for the quasi-elastic barrier distribution for the ^"Zn + 

I ^"^Pb reaction. We show that the coupled-channels calculations which include the multi-phonon ex- 

^— ( . citations in the colliding nuclei reproduce reasonably well the experimental excitation function for 

' quasi-elastic scattering at backward angles and the barrier distribution for this reaction. 
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INTRODUCTION 



! The internal structure of colliding nuclei strongly influences heavy-ion collisions at 

^ I energies around the Coulomb barrier. A well known example is a reaction of a deformed 

nucleus. In this case, the nucleus-nucleus potential depends on the orientation angle of 
the deformed nucleus with respect to the beam direction. Assuming that the orientation 
^ I angle does not change during the collision, the cross section can then be obtained by 

■ averaging the contribution from all possible angles flLQ,!!!]. In this picture, the relative 

motion between the colliding nuclei experiences many distributed potential barriers 
I depending on the orientation angle of the target nucleus, instead of a single barrier. 

' To a good approximation, the concept of barrier distribution can be extended also to 

. systems with a non-deformed target |3, where the distribution originates from 

the coupling between the relative motion and several intrinsic degrees of freedom such 
as collective inelastic excitations of the colliding nuclei and/or transfer processes. 

In Ref. [8], Rowley, Satchler, and Stelson argued that a barrier distribution can be 
directly extracted from a measured fusion excitation function afus(£'), by taking the 
second derivative of the product Eafus{E) with respect to the center-of-mass energy E, 
that is, d^{EOfufi)/dE^. This has stimulated many high precision measurements of fusion 
cross section, so that the second derivative is meaningful [4, 9]. The extracted barrier 
distributions have revealed that the concept indeed holds and the barrier distribution 
provides a nice tool to investigate the fusion dynamics of the entrance channel. It was 
also shown recently that the concept of barrier distribution still retains even in massive 
systems, such as '°^Mo -i- ^^^Mo [10]. 

A similar barrier distribution can be extracted also using the quasi-elastic scattering 



[lllill2l]. The quasi-elastic scattering is a sum of elastic, inelastic, transfer, and breakup 



processes, and is a good counterpart of heavy-ion fusion reaction [2]. A major difference 
is that the quasi-elastic scattering is related to the reflection probability of the Coulomb 
barrier, while the fusion is related to the transmission. Since the penetration and reflec- 
tion probabilities are related to each other due to the flux conservation, similar informa- 
tion can be obtained both from fusion and quasi-elastic scattering. 

In this contribution, we discuss two recent theoretical activities on heavy-ion quasi- 
elastic scattering at sub-barrier energies. We first present our recent systematic analyses 
on heavy-ion quasi-elastic scattering at deep-subbarrier energies, in aiming at extracting 



the surface diffuseness parameter of inter-nuclear potential [|13L Il4ll. We then discuss 
coupled-channels calculations for the ^'^Zn -i- ^^^Pb reaction, for which the quasi-elastic 
barrier distribution has recently been obtained experimentally 



QUASI-ELASTIC BARRIER DISTRIBUTIONS 



Before we proceed, let us first summarize the theoretical formulas for quasi-elastic 
barrier distribution. In the eigenchannel representation of the coupled-channels method, 
the fusion and quasi-elastic cross sections are given as a weighted sum of the cross 
sections for uncoupled eigenchannels [0,13, [1,0, 13]. That is. 
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where o^^J {E) and o^^' {E,d) are the fusion and the elastic cross sections for a potential 
in the eigenchannel a. Notice that the same weight factors w« appear both in Eqs. © 
and (l2l). This is a generalization of well-known orientation average formula for a system 
with deformed target. 
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where 6t is the orientation of the deformed target and represents a continuous variable 
for a in Eqs. © and ©. 

The idea of barrier distribution is led by the fact that the classical cross sections for 
fusion and quasi-elastic scattering for a single potential barrier are given by 



a',l^\E) = nRi[l-f]e{E-V,), 
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and 



of'^\E,7t) = aR{E,7t)eiVb-E) 
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respectively 11 1211 . Here, Rt, and Vt, are the barrier position and the barrier height for the 
5-wave scattering (thus the scattering angle is set to be n in Eq. ([5])), respectively, and 
Or{E, n) is the Rutherford cross section. These yield [[s, H], 
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Evidently, these functions provide information on how potential barrier heights are 
distributed, and are called fusion and quasi-elastic barrier distributions, respectively. In 
realistic situations, the quantum (tunneling) effect smears the delta function in Eqs. ^ 
and ([7]). Moreover, the effect of nuclear potential has to be taken into account in quasi- 
elastic cross sections in Eq. ([7]) il2|l . Nevertheless, from the derivation, it is apparent 
that the fusion and quasi-elastic barrier distributions behave in a similar way. This is 
demonstrated in Fig. 5 in Ref. [12] for the ^^O -i- ^^^Sm system. 

In actual experiments, it is impossible to put a detector at a scattering angle TZ. One 
can, however, scale a cross section in energy by taking into account the centrifugal 
correction. Estimating the centrifugal potential at the distance of closest approach for 
the Rutherford scattering, r^, the effective energy may be expressed as ifllll 

l^lrl ^^l + sin(0/2)- 

Therefore, one expects that the function —d/dE{oQ\/oR) evaluated at an angle d will 
correspond to the quasi-elastic barrier distribution at the effective energy given by eq. 



INTER-NUCLEUS POTENTIAL AND DEEP-SUBBARRIER 
QUASI-ELASTIC SCATTERING 

Let us now discuss the application of deep-subbarrier quasi-elastic scattering to the 
problem of surface diffuseness anomaly in heavy-ion potential For calculations 
of elastic and inelastic scattering, which are sensitive only to the surface region of the 
nuclear potential, the diffuseness parameter of around 0.63 fm has been conventionally 
employed iTylflsll. This value of surface diffuseness parameter has been well accepted, 
partly because it is consistent with a double folding potential [19]. In contrast, a recent 
systematic study has shown that experimental data for heavy-ion fusion reactions require 
a much larger value of the diffuseness parameter, rangingjpetween 0.75 and 1.5 fm, as 
long as the Woods-Saxon parameterization is employed [20]. 

Since quasi-elastic scattering and fusion are complementary to each other, it is of 
interest to investigate this problem using quasi-elastic scattering. In doing so, we are 
particularly interested in the deep sub-barrier region [[l3l [l^j. At these energies, the 
cross sections of (quasi-)elastic scattering are close to the Rutherford cross sections, 
with small deviations caused by the effect of nuclear interaction, VW. This effect can be 



taken into account by the semiclassical perturbation theory il2|,|21|], which leads to 



dOR{E,9) ka E 



where k = ^JTjlE /fi, jj. being the reduced mass, rj is the Sommerfeld parameter, and 
rc is the distance of closest approach. This formula shows that the deviation of the 
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FIGURE 1. The ratio of the quasi-elastic to the Rutherford cross sections at Oi^h = 159° for the ^^S + 
'^^Au (the upper panel) reaction and for the ^"'^S + '^^Au (the lower panel) reaction. 



elastic cross sections from the Rutherford ones is sensitive predominantly to the surface 
region of the nuclear potential, especially to the surface diffuseness parameter a. There 
is another advantage of using the deep sub-barrier data. That is, the effect of channel 
coupling on quasi-elastic scattering can be disregarded at these energies, since the 
reflection probability is almost unity irrespective of the presence of channel couplings 
[[l4]. From these considerations, it is evident that the effect of surface diffuseness 
parameter can be studied in a transparent and unambiguous way using the large-angle 
quasi-elastic scattering at deep sub-barrier energies. 

Figure 1 compares the experimental data with the calculated cross sections obtained 
with different values of the surface diffuseness parameter for the ^^S -i- ^^^Au system 
(the upper panel) and the ^"^S -i- ^^^Au system (the lower panel). In order to analyze 
the experimental data at deep sub-barrier energies, we use a one-dimensional optical 
potential with the Woods-Saxon form. Absorption following transmission through the 
barrier is simulated by an imaginary potential that is well localized inside the Coulomb 
barrier. The best fitted values for the surface diffuseness parameter are a = 0.57 ±0.04 
fm and a = 0.53 ±0.03 fm for the ^^S and ^^S -i- ^^^Au reactions, respectively. The 
cross sections obtained with these surface diffuseness parameters are denoted by the 
solid line in the figure. The dotted and the dot-dashed lines are calculated with the 
diffuseness parameter of a = 0.80 fm and a = 1.00 fm, respectively. It is clear from 
the figure that these spherical systems favor the standard value of the surface diffuseness 
parameter, around a = 0.60 fm. The calculations with the larger diffuseness parameters 



underestimate the quasi-elastic cross sections and are not consistent with the energy 
dependence of the experimental data. We obtain a similar conclusion for the 32,34g _,_ 
208p|^ and ^^O + ^ospj, systems [13]. This indicates that the double folding procedure is 
valid at least in the surface region and for spherical systems which we studied. 

For deformed systems, such as ^^O + ^^^Sm, ^^^W, on the other hand, we found that 
the surface diffuseness parameter of a = 1 . 14 ± 0.03 fm and 0.79 ±0.04 fm for the for- 
mer and for the latter, respectively, is required in order to account for the experimental 
data [|l3Lll4|] . Although these large values of surface diffuseness parameter are consis- 
tent with that extracted from fusion, the origin of the difference between the spherical 
and the deformed systems is not clear. In order to clarify the difference in the diffuse- 
ness parameter, apparently further precision measurements for large-angle quasi-elastic 
scattering at deep sub-barrier energies are urged, especially for deformed systems. 



COUPLED-CHANNELS CALCULATIONS FOR QUASI-ELASTIC 
BARRIER DISTRIBUTION FOR ^^ZN + ^o^pb REACTION 

We next discuss the barrier distribution for synthesis of superheavy elements. When 
one discusses a fusion reaction to synthesize superheavy elements, one often refers to a 
single potential such as the Bass barrier jl'^ . On the other hand, the effect of channel 
coupling is in general strong for massive systems, and thus one can expect a broad 
distribution of potential barriers. It is thus important to study how the potential barrier 
is distributed for massive systems, since it is crucial to choose the right beam energy 
in order to effectively synthesize superheavy elements. Moreover, there is no a priori 
evidence why the Bass barrier is reasonable in the superheavy region. For these reasons, 
the quasi-elastic barrier distribution measurements have been recently performed by 
Mitsuoka et al. for systems relevant to cold fusion reactions, ^^Ti, ^^Cr, ^^Fe, ^^Ni, ^^Zn 



-I- ^*^^Pb llSfl. In this section, we perform coupled-channels calculations for the '^Zn -i- 



70^ 



208p|, system. 

The calculations are done with a version [23] of the coupled-channels code CCFULL 



|24|1 . This code treats the coupling to all orders in the coupling hamiltonian and em- 
ploys the isocentrifugal approximation in order to reduce the dimension of the coupled- 
channels equations. It has been shown that the isocentrifugal approximation works well 
for quasi-elastic scattering at backward angles [12]. In the code, the regular boundary 
condition is imposed at the origin, instead of the incoming boundary condition. 

Figure 2 shows the excitation function of the quasi-elastic scattering (the upper panel) 
and the barrier distribution (the lower panel). The deep-inelastic component has been 
subtracted from the experimental data using a statistical code, as is explained in Ref. 
i\5]. The solid and dashed lines are the results of the coupled-channels and the potential 
model calculations, respectively. We use the Woods-Saxon potential with Vq = 140 MeV, 
ro = 1.186 fm, and a = 0.69 fm for the real part and Wq = 50.0 MeV, r^, = 1 .0 fm, and = 
0.1 fm for the imaginary part. In the coupled-channels calculation, we include the double 
quadrupole phonon excitations in the ^^Zn and the triple octupole phonon excitations in 
the ^'^^Pb nucleus. In addition, we include the mutual excitation channels, [1,1], [1,2], 
[2,1], and [2,2], where [np, nj] denotes the excitation channel with np phonon state in the 
projectile and nr phonon state in the target nucleus. In this way, we include 10 channels 




FIGURE 2. The ratio of the quasi-elastic to the Rutherford cross sections (the upper panel) and the 
quasi-elastic barrier distribution (the lower panel) for the ^^Zn + ^^'^Ph reaction. These are plotted as 
a function of effective energy defined by Eq. (8). The solid line is the solution of coupled-channels 
equations, which take into account the double quadrupole phonon excitations in the ^"Zn nucleus and 
the triple octupole phonon excitations in the ^"^Pb nucleus. The dashed line shows the result without the 
couplings. The experimental data are taken from Ref. [il5(l . 



(including the entrance channel, [0,0]) in the calculations. The excitation energy for the 
single phonon state and the deformation parameter are E2+ = 0.885 MeV and ^2 = 0.228 
for the projectile nucleus ^^Zn and £3- = 2.614 MeV and jh,=0.ll for the target nucleus 
^°^Pb. We use = 1.2 fm for the coupling term. 

In the code, the coupled-channels equations are solved by constructing linear inde- 
pendent solutions of the equations, where is the dimension of the coupled-channels 
equations. A linear superposition of these solutions is then taken to construct the physi- 
cal solution, which fulfills the asymptotic boundary condition for scattering. For massive 
systems, it is sometimes difficult to numerically maintain the linear independence of the 
solutions, since the wave functions scale very differently from one channel to another. 
This leads to a numerical instability of the solution of the coupled-channels equations. 
We avoid this difficulty by taking a linear superposition of the solutions at several places, 
with an interval of 1 fm up to 15 fm, so that the linear independence is recovered. See 
Ref. 1231] for details. Even though we use this prescription, we still find a small spu- 
rious oscillation in the calculated excitation function of quasi-elastic cross section due 
to the numerical inaccuracy, when the coupling is strong. We therefore average the cal- 
culated cross sections with a Gaussian weight with 0.5 MeV width. We have checked 
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FIGURE 3. The ratio of the quasi-elastic to the Rutherford cross sections (3(a) and 3(c)) and the quasi- 
elastic barrier distribution (3(b) and 3(d)) for the ^"Zn + ^*^^Pb reaction. The figs. 3(a) and 3(b) are obtained 
by including the different number of octupole phonon excitations in the target nucleus as indicated in the 
inset, together with the double quadrupole phonon excitations in the projectile nucleus. The figs. 3(c) and 
3(d) are obtained by including the different number of quadrupole phonon excitations in the projectile 
nucleus together with the double octupole phonon excitations in the target nucleus. 



that the shape of quasi-elastic barrier distribution is insensitive to the value of the width 
parameter. 

As we can see in the figure, the potential model calculation (the dashed line) sig- 
nificantly underestimate the quasi-elastic cross sections at energies above the Coulomb 
barrier. Also, the barrier distribution has a significantly narrow width, and is inconsis- 
tent with the experimental data. On the other hand, the coupled-channels calculation 
(the solid line) well reproduces the experimental data both for the excitation function 
and barrier distribution. 

Figure 3 shows the role of multi-phonon excitations. Figures 3(a) and 3(b) are ob- 
tained by varying the number of octupole phonon excitations in the target while keeping 
the double phonon excitations in the projectile nucleus. On the other hand, figs. 3(c) and 
3(d) are obtained by varying the number of quadrupole phonon excitation while keeping 
the number of octupole phonon excitations in the target nucleus to be two. These figures 
show that the double phonon excitations considerably alter the shape of barrier distribu- 
tion as compared with the barrier distribution obtained with the single phonon excitation. 
For both in the projectile and in the target nuclei, the double phonon excitation leads to 
better agreement with the experimental data, although we find that the triple phonon ex- 



citation in the target nucleus plays a much less important role. A similar conclusion has 
been obtained also in Ref.[10], where the role of multi-phonon excitations was discussed 
for the ^^*^Mo+^''^Mo fusion reaction at energies around the Coulomb barrier. 

SUMMARY 

We have discussed two recent developments in heavy-ion quasi-elastic scattering at ener- 
gies around the Coulomb barrier. We first discussed the surface property of internucleus 
potential. We have argued that the quasi-elastic scattering at deep subbarrier energies 
offer a clear and almost model independent way to determine the surface diffuseness 
parameter, that is, the slope of asymptotic exponential tail of the potential. The value of 
diffuseness parameter extracted from the 3234g _^_ 197^^ reactions is around 0.55 fm, and 
is consistent with the double folding potential. On the other hand, the surface diffuseness 
parameter extracted from systems with a deformed target, that is, ^^O -i- ^^"^Sm, ^^^W 
was found to be much larger (1.14 fm for the former and 0.79 fm for the latter systems). 
Further investigations will be required in order to clarify the system dependence of the 
surface diffuseness parameter. In the second part, we performed the coupled-channels 
analyses for a cold fusion reaction ^^Zn -i- ^*^^Pb, where the quasi-elastic barrier distribu- 
tion was recently obtained by Mitsuoka et al. . Including the double quadrupole phonon 
excitations in the projectile nucleus '^^Zn and the triple octupole phonon excitations in 
the target nucleus ^''^Pb in the coupled-channels calculation, we could reproduce rea- 
sonably well both the excitation function of quasi-elastic cross section and the shape of 
quasi-elastic barrier distribution. This indicates that the coupled-channels approach still 
works for the approaching phase of the reaction even in massive systems, where many 
degrees of freedom may be involved in the reaction [10]. It also suggests that the deep- 
inelastic collision can be regarded as a post-barrier phenomena, since the experimental 
quasi-elastic cross sections have been obtained by subtracting the deep-inelastic compo- 
nents. We also discussed the role of multi-phonon excitations, and showed that they play 
an important role in this system. The coupled-channels analyses for other cold fusion 
reactions, ^^Ti, ^"^Cr, ^^Fe, ^"^Ni -i- ^'^^Pb are now in progress, and we will report on them 
in a separate publication. 
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